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Abstract 

We examine the possibility that the form of the lepton mixing matrix can be 
determined by extremising the Jarlskog flavour invariants associated, eg. with 
the commutator (C) of the lepton mass matrices. Introducing a strictly covariant 
approach, keeping masses fixed and extremising the determinant (Tr C^/3) leads 
to maximal CP violation, while extremising the sum of the 2x2 principal minors 
(— Tr C'^/2), leads to a non-trivial mixing with zero CP violation. Extremising, 
by way of example, a general linear combination of two CP-symmetric invariants 
together, we show that our procedures can lead to acceptable mixings and to 



non-trivial predictions, eg. \Uez\ — \/2/3 y Amf2/Am23(l — m^/nir)'^ ~ 0.07. 
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1. Introduction 

It is now almost twenty years since Jarlskog [1] first alerted us to the inherent U(3) 
invariance representing the freedom in any weak basis [2] to transform fermion mass 
matrices in the generation space, whilst keeping the charged-current weak-interaction 
diagonal and universal, ie. whilst always staying in a weak basis. While all observables 
(eg. masses and mixing angles) must of course be Jarlskog- invariant, we focus here on 
fiavour-symmetric invariants, ie. those invariants without fiavour indices of any kind 
(ie. involving no flavour projection operators in their definition). 

While it has long been clear [3] that any fundamental laws underlying the masses 
and mixings should be Jarlskog-covariant (as defined below), this principle may not 
always be easy to respect in practice. In this paper we will show that cxtremisation 
of flavour-symmetric Jarlskog invariants [4] [5] leads naturally to Jarlskog-covariant 
constraints, which (as we shall see) can even have viable phenomenology. 

The archetypal example of a flavour-symmetric Jarlskog invariant is the famous 
determinant of the commutator [6], which for the charged- lept on (L) and neutrino 
(iV) mass matrices (taken here to be hermitian) may be written: 

DetC = TrCV3 = Tri[L,N]y3 = -2 Det diag(A,) Det diag(A^) J (1) 
Af ^ {mi^-nir, rrir -nie, me-m^), ^ {m2 - m^, ms - mi, mi - 1112) (2) 

where C := —i[L,N], with me, m^, mr the charged-lepton masses and mi, m2, m^ 
the neutrino masses, and with the invariant measure oi CP violation [6] given by: 

J = C12 S12 C23 S23 Si3 Ss = {I - 5^2)^ S12 (1 - 533)^ "523(1 " "Sl3 83. (3) 

Extremising Eq. 3, with respect to the standard PDG [7] variables S12, S23, ^13 and S is 
known [8] to lead to S12 = S23 = l/\/2, S13 = l/\/3 and S — n/2, corresponding 

to maximal CP violation, ie. to trimaximal mixing [9]. 

It should be remarked that the assumption that the mass matrices are hermitian 
(see above) may in fact be realised in two distinct ways: one may simply imagine 
applying a suitable transformation to the right-handed flelds (which are anyway inert 
to the charged-current weak-interaction) bringing L and N into hermitian form, or one 
may instead everywhere re-interpret L and N as hermitian squares of masss matrices 
L — > LL^, N — > NN'^, operating between left-handed fields. In the latter case, masses 
need to be replaced by masses-squared throughout {mi/i, — > mf^^, eg. in Eqs. 1-2). 
In this paper, in order to circumvent this as yet unresolved ambiguity, we shall give 
results, where appropriate, for both the linear and quadratic cases. 

Of course we now know that lepton mixing is not trimaximal, it being actually 
much closer to the so-called tri-bimaximal form [10] [11]. The question then arises, as 
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to whether extremising some other invariant (or combination of invariants) aside from 
TrC^, might perhaps yield a more reahstic mixing prediction. 

In this paper we seek to kindle interest in these kinds of questions, and we start 
out (Section 2) by re-deriving the results for Tr above, introducing and bringing to 
bear a strictly covariant (basis-independent) approach. We focus here on the leptonic 
case: in the context of extremisation, quasi-maximal mixing in the lepton sector looks 
a priori more comprehensible than anything involving the quarks (which are anyway 
subject to larger radiative corrections, eg. dependent on the top mass). 

We proceed (Section 3) with a similar covariant extremisation of the only other 
independent invariant depending on C alone, ie. TrC^. The mixing prediction is now 
CP-conserving and, whilst not phenomenologically viable as it stands, is non-trivial, 
with even a suggestive resemblance to some previously proposed Icpton-mixing ansatze. 

Finally, taking a particular linear combination of two CP-symmetric invariants as 
an illustrative example, we show (Section 4) how covariant extremisation of slightly 
more general combinations of invariants can readily lead to rcalisic mixing predictions. 

Our "Perspective" section (Section 5) simply consolidates what we have learned, 
also briefly pointing forward to what one might hope to learn in future studies, eg. in 
terms of incorporating CP violation, dropping mass constraints etc. 

2. Extremising the Jarlskog Determinant, Tr 

One might argue that the PDG variables are arbitrary, eg. it matters (Eq. 3) that we 
chose to extremise with respect to 5 and not ss — sin 5. We choose to extremise with 
respect to the Yukawa couplings themselves (ie. with respect to the mass matrices) 
using some established results [13] on differentiation of matrix traces with respect to 
matrix variables (see Appendix A). In terms of matrix derivatives Sl '■— d/dL etc., 
the extremisation conditions for, eg. Tr C^ may then be written (see Appendix A) in 
a manifestly Jarlskog-covariant (ie. basis-independent) way: 



in that both sides of the equation are evidently form-invariant under an arbitrary 
C/(3) transformation L C/(3)LC/(3)^ N U{?>)NU{?>)\ ie. we would obtain the 
same equations performing the extremisation in any basis. When the masses are to 
be held fixed, the zeros on the RHS of Eqs. 4-5 must be replaced by arbitrary matrix 
polynomials in and A^-^ respectively, representing the dependence on unknown 
Lagrange multipliers (see below), with the manifest covariance clearly maintained. 



Ol Tr CVS 
dN Tr CV3 



-i[N,CY = Q 
+i[L,C2]^ = 0, 



(4) 
(5) 
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The covariance property allows us to solve our equations in any basis, so we choose 
a basis convenient to us where the charged-lepton mass matrix is diagonal and where 
all imaginary parts of the off-diagonal elements of the neutrino mass matrix are equal 
(the 'epsilon' basis [12]). The neutrino mass matrix may then be written: 



(6) 



where the seven variables a, fe, c, x, y, z and d determine the three neutrino masses 
and the four mixing parameters (the charged-lepton masses being directly the diagonal 
elements of L). In terms of these variables the Jarlskog determinant is given by: 

Det C = Tr = 2d{me — m^)(m^ — mT){m.j- — me){(f — xy — yz — zx) (7) 

so that CP violation vanishes in the case d = and also in the case d'^ = xy-\-yz + zx. 

The off-diagonal elements of —i[N, C^] must be set to zero for an extremum (Eq. 4), 
leading to the following cyclically symmetric constraints, from the real (Re) parts: 



-d {me - m^) {rrir 


-nie) ({y-z) - 


(b- 


■c)) {y + z) = 







Re — d {m^ — rrir) {f^e 


-rrifj) (y{z-x) - 


{c- 


a)) {z + x) = 







-d {rrir - rrie) {m^ 


-rrir) {{x-y) - 


{a — 


■b)) {x + y) = 







and from the imaginary (Im) 


1 parts: 










-{rUe - m^) {rrir - 


■me) ((6-c)(d2- 


yz) 


- {y'-z')x) 







Im — {rrin — rrir) {'^e — 


m^) ( (c - a) {d^ - 


zx) 


- {z'-x')y) 







-{rrir - rrie) (m^ - 


rrir) {{a-b){d^- 


xy) 


- {x' - y') z) 




0. 



(8) 



(9) 



It should be remarked that the cyclic symmetry (e— >t, a— >6— >c, x— >y— ^z) 
will be seen to be a useful and important generic feature of our approach, resulting 
from starting with flavour-symmetric invariants and working in the epsilon basis. 

A non-trivial solution to Eq. 8, ie. that corresponding to the real parts above 
(putting aside the d — possibility for one moment), may be written: 

{x — y) = {a — b) a — a + X 

{y-z) = {b-c) i.e. b^(7 + y (10) 

{z — x) — {c — a) c — (7 -\- z 

where a is an undetermined overall constant offset. Eq. 10 is precisely the 'S3 invariant 
constraint' [14] (or 'magic square constraint' [12]) whereby the neutrino mass matrix 
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a = b 


and 


X = y 


b^c 


and 




c — a 


and 


Z — X 



is determined to be an 'S3 group matrix' [14], ie. a matrix having all row/column sums 
equal and hence (at least) one trimaximal eigenvector. 

Turning to Eq. 9 (corresponding to the imaginary parts) a non-trivial solution is: 



(11) 



whereby the mass matrix is circulant [4] and the mixing takes the trimaximal form [9] 
(Eq. 10 being also satisfied). More trivially, coming to the d — case (above), we have 
that twofold-maximal mixing, eg. b — c, y — z — also satisfies both sets of equations. 

Of course, these solutions do not solve all the extremisation equations, Eqs. 4-5. It 
transpires however, that in the case that the masses are held fixed, all the remaining 
equations serve only to determine suitable Lagrange multipliers Xu, Ajvi {i — 0, 1, 2): 

(dLTrCy3f = -t[N,C^] : 
{dNTrCysf = +i[L,C^] : 



Alo + Ali L + Al2 L 



(12) 
(13) 



with the apparently excess constraints turning out to be redundant (and furthermore 
with this circumstance occuring for all the cxtrcmisations studied in this paper). 

In the epsilon basis, L is diagonal and any polynomial in L is diagonal also, whereby 
Eqs. 8-9 (and their solutions, eg. Eqs. 10-11) remain valid, even with the Lagrange 
multipliers (Eq. 12). In addition, the on-diagonal elements of +i[L,C^] vanish (as 
in Eq. 5) due to L being diagonal, and in the case of non-zero Lagrange multipliers 
(Eq. 13) the consistency requirement on the coefficients of the Xm takes the form: 



a 
b 
c 



a? + y'^ + 
x"^ + b'^ + z^ 
+ + ^ 



0. 



(14) 



All our solutions above do indeed turn out to satisfy this determinant condition. 
(In fact any 5*3 group matrix, Eq. 10, automatically satisfies this condition, Eq. 14). 

Covariance implies that the XujNi themselves will be Jarlskog scalars, expressible 
in terms of, eg. traces of powers of mass matrices: 



= TrL = 
= TrL^ 



rrie + rrin + 



m, 



ml + + ml 



TiN = 
TrAT^ 



mi + 7722 + 



mi 



mo 



mq 



(15) 



ml + ml + ml 



or traces of commutators of mass matrices etc. Note that the definitions (Eq. 15) are 
our 'constraint equations', eg. TrL^ = L2 for fixed L2 etc., which when differentiated, 
Ol (Tr - L2) = 2L^ etc., lead to the polynomial forms on the RHS of Eqs. 12-13. 
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The on-diagonal elements of Eq. 12 (which are purely real) are given by: 



2d (m^ + mr- 2me) {m^ - rrir) {d - x - y - z ) ^ Xlq + XLiUie + XL2m^ 

2d {rrir + me - 2m^) {rrir - me) {d^ - - - z^) = A^o + A^im^ + Ai2m^ (16) 

2d (me + m^ - 2m^) (me - m^) (d^ - x^ - y"^ - z^) = A^o + A^im^ + Ai2m^. 

Solving for the Xu gives the charged- Icpton discriminant L\ in the denominator of the 
solution (from the determinant-of-coefficients, Eq. 16, RHS), where La is given by: 

La := ^Li/2 + 3LfL2/2 + 6L1L2L3 - 7LfLi/2 - 3Li - AL\L^/^ - L?/6 

= {rUe- m^){m^- mr){mr - me). (17) 

After some work to cast the numerators also into invariant form, we find that: 



-Lf /2 + 3L?L2 - 7LiL2/2 - 2L2L3 + 3L2L3 
3Lt/2 - 7L2L2 + 3L2/2 + 6L1L3 



A^o = ^^1^3 ^^^2^3 ^^3 (18) 



Ail = " -i"- — "^/- ■ o (19) 

-9L3 + 9L1L2 - 2L? , , 

Al2 = ^-TiC\ (20) 

Of course, these expressions can also be obtained in a fully covariant way simply by 
multiplying Eq. 12 by successive powers of L (/ = L", L = L^, L^) and taking traces 
(powers higher than L^ must always be reduced using the characteristic equation). 
Analogous L ^ N expressions for the X^i solve all the remaining equations. 

Recognising the coefficients of TrC^ on the RHS above as dL/^iL /^N js) / {L /^N /:^) , 
we see that we have, in effect, extremised the Jarlskog invariant J = TrC^/(LAiVA), 
recovering the known result [8], and validating our procedure. Finally, we note that 
the above analysis is essentially unchanged, taking L and N as the hermitian squares 
of mass matrices (see Section 1) substituting masses-squared for masses throughout. 

3. The Sum of the 2x2 Principal Minors, -TrC^/2 

The other independent invariant function of C may be taken to be the sum of the 2x2 
principal minors, expressible in terms of the X-matrix [Ki^ :— —Ki,\Yi" \lb\) as follows: 

Qn = -Tr C /2 = Tr [L, Nf /2 = Af diag(Ai) K diag(A,) A,. (21) 

The matrix is a key oscillation observable [15], carrying equivalent information to 
the moduli-squared of the mixing elements [16]. The matrix comprises the real parts 
of the plaquette products [17]: 11;,^ := —Ki^, -\- iJ (indices to be interpreted mod 3) 
and may be viewed as the CP-conserving analogue of the Jarlskog invariant J [6]. 
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Prom Eq. 21, taking for simplicity initially the extreme hierarchical approximation 
rrie, m^, mi, — > (ie. the 2x2 mixing limit), we have: 

= mlml c^gc^a (si3 + c^gS^g) + 0{mlml) + ... (22) 

Extremising this expression with respect to S23 and S13 gives S23 = S13 = or 
"S23 = — 1 or C23C13 = l/\/2 (where the latter looks at first sight very promis- 
ing phenomenologically, see below). Extremising rather with respect to ^23 and ^13 
gives in addition S23 = 1 or S13 = 1 as solutions, for arbitrary S13 and S23 respectively. 
Since mixing matrices which are just permutation matrices give no mixing as such, we 
have that 2x2 maximal mixing is the only non-zero mixing solution in the 2x2 case. 

Just as for Tr above, however, we will extremise exactly here with respect to 
the mass matrices (see again Appendix A): 

{-dLTTCy2f = +[N,[L,N]] = 0, Xl, + XliL + Xl2L^ (23) 
{-dNTrCy2f = -[L, [L,N]] = 0, Xm + XmN + XmN^ (24) 

where the matrix polynomials replace the RHS zeros, in the case of mass constraints. 
In the epsilon basis the off-diagonal elements of Eq. 23 are, for the real parts: 

(m^ + rrir — 2me) {(f — y z) + (m^ — rrir) {b — c) x = 
Re {rrir + me — 2mfj_) {(f — z x) + {mr — me) {c — a) y = (25) 
{me + m^ — 2m^) {d^ — xy) + {me — m^) {a — b) z — 

and for the imaginary parts: 

d (m^ -|- mr — 2me) {y + z) + d {m^ — mr) {b — c) — 
Im d {mr + me — 2m^) {z + x) + d {mr — me) {c — a) — (26) 
d {me + mn — 2mr) {x + y) + d {me — m^) {a — b) — 0. 

We observe that Eqs. 25-26 imply that either d = or = xy + yz + zx, and so from 
Eq. 7 we see immediately that there can be no CP violation in this case. 

It turns out that the d'^ = xy+yz+zx solution is equivalent to the d = solution, so 
that we need to consider only the d—0 case (note that in the case — xy+yz + zx the 
epsilon basis is not unique, and the imaginary part can be rephased to zero). Setting 
d — then, solves Eq. 26, corresponding to the imaginary parts. Clearly, setting any 
pair of X, y, z to zero and the corresponding pair of a, b, c equal, solves also Eq. 25 
for the real parts, so that, just as for Tr above, 2x2 maximal mixing in any sector 
gives an exact extremum of Tr C^, independent of the neutrino mass spectrum. 
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A somewhat more interesting non-trivial solution to Eq. 25 (still with d — from 
Eq. 26) clearly reflects the inherent cyclic symmetry: 



X — ±y (a 


-b){c 


-a)E 




-c){a- 


-h)M 




-a){b- 


-c)T 



(me 


- m^){mr 




) 


(me + - 


- 2m,-) (m^ 


+ "^e 


- 2m/,) 


{m^ 


- m^)(me 




) 


(m^ + mr- 


- 2me)(me 


+ 


- 2m^) 


{rrir - 


- me) (m^ 






{rriT + me — 


2m^)("^/. 


+ rrir ■ 


- 2me) ■ 



E 

M = ^""^ """O'""" ""^^ ^ (27) 

r = 

With the charged- Icpton masses positive, we have > > M, T, whereby cither 
h < a < c or c<a<6. We may take all upper (positive) signs in the case 
6 < a < c, corresponding to an inverted hierarchy, and all lower (negative) signs 
in the case c < a < 6, corresponding to a normal hierarchy (as considered below) 
with other sign combinations being just trivial rephasings of these two possibilities. 
For these solutions the determinant condition is again automatically satisfied, since 
£'-|-M-|-T'-|-l = appears as a factor in the determinant Eq. 14. 

As in the case of Tr above, the remaining constraints (from Eqs. 23-24) will serve 
only to fix the Lagrange multipliers. While we again have the option to determine the 
Lagrange multiplers in a general and fully covariant way, it will prove useful as before 
to consider eg. the on-diagonal Ol equations (which are again purely real) explicitly: 

2{mr-me)y'^ - 2(me-m^)2;^ = A^o + A^ime + AL2"^e 

2{me- m^)z^ - 2{mij, - mr)x'^ = Alo + ^^Lirrifj, + Al2"^^ (28) 

2(m^-mT-)a;^ - 2(m^-me)y^ = A^o + A^im^ + \L2ml. 

Since Lagrange multipliers need in general only be evaluated locally at the extremum, 
we now have the possibihty to achieve some simphfication by substituting (on the LHS 
of Eq. 28) the specific solutions Eq. 27, before solving for the Xu and again casting 
the resulting expressions (curly brackets denote anticommutators) in invariant form: 



(Tr{L2,7V}Li - Tr {L, A^} L2)(3Tr {L , TV} - 2L^Ni) 
2(9L3 - 9L1L2 + 2L\] 



AiO = nr r , or3N \^^) 



_ -(3Tr{L^A^}-2L2A^l)(3Tr{L,A^} - 2LiAr,) 
-^^^ " 2(9L3 - 9L,L2 + 2L3) ^"^"^ 

^ (3Tr{L,iV} - 2L,N,f 

2(9L3 - 9L1L2 + 2L?) ■ ^ 

(Similar expressions also result for L ^ A^, as before). While this means that our 
Lagrange multipliers are now solution-specific (ie. they fail for the mass-independent 
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2x2 maximal-mixing solution above), one at least has that Eqs. 23-24 and Eqs. 29- 
31 (with their L N counterparts) together constitute a covariant statement of 
specifically the non-trivial (cyclically symmetric) solution given by Eq. 27. 

Prom Eq. 27 the mixing is uniquely predicted, in essence by fitting a, b and c (Eq. 6) 
to give the correct neutrino mass spectrum. Of course absolute neutrino masses (as 
distinct from mass-squared differences) are not yet known directly, but we may still 
obtain a specific prediction if we are prepared to assume a normal (ie. non-inverted) 
'classic' fermion mass hierarchy mi << m2 << ms (ie. with no offset, toi ~ 0). 

Wc then find that an acceptable neutrino mass-squared difference hierarchy {hi := 
Am^2/'^^23 — "'l/'^'l — 0.03) results, simply by setting the ratio {b — a)/{a — c) ~ 1, 
corresponding, it turns out, to the Ue being exactly trimaximally mixed (this ratio is 
in fact the only operative parameter here, since rescaling and off-setting of a, b and c 
clearly does not inffuence the mixing matrix). The resulting |f/p-matrix (ie. the matrix 
of moduli-squared of the resulting mixing-matrix elements) is given numerically by: 



















J^2 




e 


f 


.33333 


.33333 


.33333 


\ 


e 


( 1/3 


1/3 


l/3\ 






.17079 


.16257 


.66663 






1/6 


1/6 


2/3 


T 




.49587 


.50409 


.00003 


) 


T 


I 1/2 


1/2 


/ 



While being closely a transpose/permutation of the tri-bimaximal form, and bearing 
a clear relation to some previously proposed ansatze (see in particular Acker et al. 
[18] and Karl et al. [19] as well as the Pritzsch-Xing ansatz [20]), this mixing (Eq. 32) 
clearly cannot at this point be considered acceptable phenomenologically. 

We find that the situation is not improved substituting masses-squared in place of 
masses throughout. Indeed the large discrepancy in \Ur3\ is considerably worsened, 
with |[/r3p ~ 3 X 10^^ (maintaining the mass hierarchy unchanged requires setting 
{b — a)/{a — c) ~ 0.36). Having also explored less minimalist assumptions regarding 
the neutrino mass-spectrum, ie. inverted mass hierarchy, large mass offsets etc. we can 
report no phenomenlogically acceptable solutions to the problem of extremising Tr C^. 
For mi << 1712 << fns, we see that all our exact solutions are consistent with the 
(approximate) constraints derived differentiating Eq. 22 with respect to PDG angles 
(although we learn that the initially promising Urs = C23C13 = 1/ condition in fact 
corresponds only to phenomenologically uninteresting 2x2 maximal mixing solutions). 

4. Extremising More General Mixing Invariants 

More general commutator invariants are readily constructed, involving higher powers 
of mass matrices [21]. In this section we shall extremise a very simple composite 
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function, comprising an arbitrary linear combination of two such invariants: 

A := Qn + gQsi (33) 

with q a scalar constant of suitable dimensionality (assumed Jarlskog-invariant) and: 

Qn := TV [L, N] [L, N] / 2 Qsi TV [L, TV] [L\ N] / 2, (34) 

our Qmn notation being essentially self-explanatory ^ at this point. 

Of course, Qn has already been individually extremised in Section 3 (Eqs. 23-24), 
so we only need in addition here the relevant derivatives of Q21 (see Appendix A): 

{OQ^ildLf = [iV, [L^ N]]/2 + {L, [AT, [L, Ar]]}/2 (38) 
{dQ^i/dNf = -[i^',[i^,Ar]] (39) 

where again curly brackets signify the anti-commutator. 

The off-diagonal terms derived from (521, analogously to Eq. 25, are (Re parts): 

(m^ -|- -|- m^m,- — rrirme — meUin — ml) {d^ — y z) + (m^ — m^) (b — c) x 

{ml + ml + m-j-me — mem^ — m^m-r — m^) {(f — zx) -\- (m^ — ml) {c — a)y (40) 

{ml + m^^ + mjififj, — m^mr — m^me — m^) {(f — xy) + {ml — m^) {a — b) z 

and (Im parts): 

d (m^ + ml + mfjjnr — m^me — mem^ — ml) {y + z) + d (m^ - ml) {b — c) 

d {ml + ml + mT-me — mem^ — m^mr — m^) {z + x) + d {ml — ml) (c — a) (41) 

d {ml + -|- mgm^ — mfjrir — m^-me — ml) {x + y) + d {ml — m^) (a — b) 

and the combined constraint equations are obtained simply by adding these terms 
(weighted by q) into Eqs. 25-26. 

^We are considering here flavour-symmetric quadratic commutator invariants (vanishing in the 
case of zero mixing) which may be usefuhy arranged as a 3 x 3 matrix: 

/ Tr[L,7V]2 Tt[L,N][L,N^] Tt[L,N^]^ \ 

Q = - { Tt[L,N][L^,N] Tr[L,N][L^N^] Tt [L , N^][L\ N^] . (35) 
\ Tr[L2,iV]2 Tt[L'^,N][L'^,N'^] Ti[L'^,N'^f J 

As usual, powers of mass matrices higher than L^, are not considered, by virtue of the char- 
acteristic equation (and Tr [L'^,N] [L,N-^] = Tr [L,N] [L'^,N^], so there are indeed only nine such 
invariants). Generalising Eq. 21, the Qmn {m,n = 1,2,3) are the double moments of the K-matrix: 

Qmn = Af diag(Ai) (diag S;)""' K (diag S,)""! diag(A,) A, (36) 
T,i = {mfj_ + mr, rur +me, me + nifj,) T,,y = {m2 + m^, ms + mi, mi + m2). (37) 

For known masses, the Q-matrix and the X-matrix are therefore equivalent, and both are equivalent 
to the |J7;!yp-matrix (see [15] [16]), with only four elements functionally independent in each case. 
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As in Section 3 we may again set d = 0, and the solutions to these combined 
equations are then readily written down from Eqs. 27, by modifying the expressions 
for E, M and T in Eq. 27, making the (cyclically-symmetric) substitutions: 



(me - m^) 
(m^ - rur) 
{rrir - me) 
{nie + - 2mr) 
{m^ + mr- 2me) 
{rrir + me - 2m^) 



(42) 



me - m^) + q{ml - m^) 
rrif, - rUr) + q{ml - ml) 
mr — me) + g(m^ - m^) 

^e + m^ — 2m,-) + q{ml + m^ + mem^ — mrLi) 

m^ + mr — 2me) + q{m'i + m^ + m^mr — meLi) (43) 



m. 



2m^) + q{ml + ml + mrme - m^Li). 







The determinant condition Eq. 14 remains valid since the equation E -\- M -\-T 
(see text just after Eq. 27) itself remains valid under these substitutions. 

Fitting to the 'mass' spectrum as in Section 3, we now have (replacing the unique 
mixing prediction Eq. 32) a trajectory of mixings, each mixing given by a different 
value of the parameter q. To locate the phenomenologically viable solution along this 
trajectory (supposing even that there is a viable solution), we will simply impose the 
"5*3 constraint" (Eq. 10), which although not forced on us here by the extremisation, 
is well-known [12] [14] to be consistent with the phenomenology. We shall see below 
how the 5'3-constraint (by its nature) can be imposed in a fully covariant way. 

In practice, we shall substitute the S?> constraint Eq. 10 into the solution Eq. 27 
(as modified by the substitutions Eqs. 42-43), thereby fixing all the parameters in the 
neutrino mass matrix, as functions of g, as follows: 



X -\- a 



y + o- 



— z + a 



X 



z — 



(m^ - mr)^(l + q{m^ + mr)) 

^ 

(m-r — me)^(l -l- q{mr + me)) 

1^ 

(me - m^)2(l + q{me + m^)) 



(44) 



where the mixing clearly cannot depend on the scale factor k or the offset a. 

Assuming a normal classic neutrino mass hierarchy as before, the operative param- 
eter is now q, fixing both the hierarchy ratio :— Am^2/Am23 ~ m^/m^^ and the 
mixing simultaneously. We see immediately from Eq. 44 that, eg. as g — > — l/(m^-|-m-r) 
we have a pole where x — > — oo while the hierarchy factor /i,^ — > 0, with the mixing 
approaching the exact tri-bimaximal form [10] [11] in the limit. 

Parametrising deviations from this pole in the form g = — (1 -|- e)/(m^ -|- m^), we 
find: e ~ 4/3 /if/i^(l + /i^/3 + S/i^/g + . . .). For hi ~ 0.03 (with hi := m^/m^ ~ 0.06) 
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we have e ~ (0.03)^, and taking k ~ 0.38 meV^/^GeV^/^ and cr ~ 25 meV, gives a 
neutrino mass spectrum (mi, m2, rris) ~ (0, 8.7, 50) meV, and a mixing matrix: 





^1 




^3 






^1 


1^2 


^3 




0.662 


0.333 


0.005 \ 




f 


2/3 


1/3 







0.219 


0.333 


0.448 






1/6 


1/3 


1/2 




0.120 


0.333 


0.547 j 




[ 


1/6 


1/3 


1/2 



(45) 



The mixing Eq. 45, being very close to the tribimaximal form, is certainly acceptable 
phenomenologically. There is the non- trivial prediction \Ue3\ - K{l-hif ~ 0.07 

(cf. Ref. [12]) and violations of — r symmetry [22] also, eg. 623 — 7r/4 — |?7e3|/\/2- 
Related poles in y and for which the mixing matrix would have rows correspondingly 
permuted with respect to Eq. 45, would seem not to be phenomenologically relevant. 

Repeating the analysis with masses-squared, in effect takes us closer to the pole 
(now e ~ 5 X 10~^) yielding a mixing correspondingly closer to the tri-bimaximal form: 



{\Uw?) 





^1 


1^2 


1^3 




^1 


^2 


^3 




0.6665 


0.3333 


0.0002 ^ 




( 2/3 


1/3 







0.1768 


0.3333 


0.4898 




1/6 


1/3 


1/2 




0.1567 


0.3333 


0.5100 ) 




I 1/6 


1/3 


1/2 



. (46) 



The corresponding prediction for |f/e3| (now If/esI ~ V2/3 hi ~ 0.014) is clearly of less 
immediate phenomenological interest than that in the unsquared case above. 

Lagrange multipliers for Eq. 33 are found in the usual way. For the di equations, 
solving the analogue of Eq. 28 (after the substitutions Eq. 42), we find: 

' 2iq^Lpe + 2qLp2Lps + L%,y ^ 

' 2{q'Lp, + 2qLp2Lp, + Lj,,)^ ^ 

' 2(g2Lp6 + 2gLp2Lp3 + Ll^Y ^ ' 

where we have explicitly assumed the '5*3 constrained' solution Eq. 44. Fully invariant 
and flavour-symmetric, the expressions for A^oo; Alqi • • • Alio etc., as functions of the 
Lj, i = 1 . . . 3, are given in Appendix B. 
Similarly, for the equations we find: 

X AT + 2gLp2Lp3 + g^Lp6 , . 

Aato = -N2 --— —— — ^ (50) 

6iVi(Lp2 + qLps) 

Xm = (51) 

L'^P2 + 2qLp2Lp3 + q'^Lpe 

2Ni{Lp2 + qLp3) 
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with no clear relation to the Sl case, since our "action" (Eq. 33) is not L N sym- 
metric. The supplementary polynomials Lp2, Lp^, LpQ are defined in Appendix B. 

It should be emphasised that (up to the apparently spontaneous choice of pole) the 
extremisation conditions on the composite action Eq. 33, together with the particu- 
lar Lagrange multipliers defined by Eqs. 47-52 and Appendix B, constitute an entirely 
covariant and fiavour-symmetric specification of the realistic, non-trivial mixing Eq. 45. 

5. Perspective 

Motivated by the notion (Section 1) that the fundamental laws underlying the fermion 
masses and mixings might come from a variational principle, applied to some as yet 
unspecified function of fiavour-symmetric Jarlskog invariants, we have presented a 
manifestly covariant machinary for extremising such invariants as functions of the 
mass matrices. (Building on the analogy [4] with the field-strength tensor in terms 
of covariant derivatives, cf. F^j^ = — i[D^,P^], we see our extremisation equations, eg. 
Eqs. 23-24 as analogous to the Yang-Mills [23] equations, cf. \D^, F^'^] = 0, which are 
themselves derivable from a quadratic Lagrangian C — — Tr F^/2, whereby we see the 
Yukawa couplings as "dynamical variables" , somewhat analogous to the gauge fields) . 

Having tested-out our methods on the familiar case of Tr (Section 2) we have 
been able to show (Section 3) that extremising Tr does not lead to physically 
realistic lepton mixings. Focussing on quadratic "actions", we see Eq. 33 (Section 4) 
as just an example of an (approximate) "effective" action, whose main particular merit 
is that wc know how to solve its resulting extremisation equations analytically, and 
that it leads to at least one mixing (Eq. 45) which is compatablc with observation. 
Clearly the complete absence of CP violation is a significant deficiency of Eq. 33 
(given the sucesss of the standard explanation of CP violation in the quark sector) 
and incorporating CP violation is an obvious challenge for the future. 

Most of all, we would like to think that this paper will stimulate the search for 
a natural and beautiful action function, determining not only the mixing but also 
the (relative) mass values mi : m2 and mg : : m^, obviating the need for 

those ugly Lagrange multipliers (and perhaps also shedding some light on the quark 
case at the same time). Regarding masses, it is clear that our procedures apply 
directly as they stand: one will only have to notice that for some 'perfect' action, 
not only is the mixing correctly predicted, but the Lagrange multiplier functions, on 
the RHS of the extremisation equations, vanish numerically for the correct masses, 
effectively as in a free extremisation. (Empirical relations like the Koide relation: 
512L3L1 - QALl - 656L2L? + 207Lf = [24] could be relevant here, cf. Appendix B). 
In removing a large part of the arbitrariness inevitably associated with the fiavour 
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sector as we know it today, such an "action" (if it exists) would surely be welcomed 
as an economical adjunct to the present Standard Model of Particle Physics. 
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Appendix A 

We are making use of a long-established result from matrix calculus [13], whereby, 
in terms of the matrix derivative dx — d/dX, for some constant matrix A, we have: 

dx Tr XA = Tr {dxX)A = A^ (53) 

where the superscript T denotes the matrix transpose. Note that dxX is a matrix of 
matrices (each sub-matrix with a single unit entry), and that in particular dxX ^ /, 
where / is the identity matrix (while (9x Tr X = Tr (9x ^ = /) • 

We first consider Tr = Tr i[L, N]^ (where C := —i[L, N], see Sections 1-2): 

^^TrC^ = Tr a^C^ (54) 

= 3Tr {dLC)C^ (55) 

= -3iTr [(dLL),N]C^ (56) 

= -3iTr {dLL)[N,C^] (57) 

= -3i[N,CY (58) 

making use of the cyclic property of the trace, Eq. 54-55 and Eq. 56-57, and the 
matrix-calculus theorem (Eq. 53) for the final step Eq. 57-58. 

The case of Tr = — Tr [L,N]'^ (see Section 3) follows entirely analogously: 

diTrC^ = -2i[N,Cf (59) 
= -2[N,[L,N]f (60) 

where in the last line we have removed the factor i using C = —i[L,N]. 

Finally, we consider the case corresponding to Q21 — Tr [L, N] [L^, N]/2 (Section 4): 

dL Tt[L,N][L^,N] = Tr [[{dLL),N][L',N] + [L, N]{[{dLL)L, N] + [L(a^L), iV])) 

(61) 

Taking each of the above three terms one at a time, we have: 

Tr [{dLL),N][L^,N] = Tr {dLL)[N,[L\ N]] (62) 

Tr [L,N][{dLL)L,N] = Tr {dLL){LNLN - L'^N'^ - LN^L + LNLN) 

= Tr {dLL){L[N,[L,N]]) (63) 

Tr [L,N][L{dLL),N] = Tr {dLL){NLNL - LN^L - N'^L'^ + NLNL) 

= TV {dLL){[N,[L,N]]L) (64) 

wherby adding the above equations: 

dLTr[L,N][L',N] = Tr {dLL){[N,[L' , N]] + {L,[N,[L, N]]}) 

= [N,[L\N]f + {L,[N,[L,N]]f (65) 

where the curly brackets denote the anti-commutator (see main text Eq. 38). 



15 



Appendix B 

The following completes the expressions for the Lagrange multiplers of Section 4: 



^Loo — Lp2{12L^Li — 6L3L2L1 — I8L3L2 

+ALsLl/3 + 19LlLi - llL^Lf + SLaL^ - L^/3) (66) 
Alio = Lp2(-36L^ + 36L3L2LI - 4L3L? - 3L^ + 3L^L? - SLaL^ + L?) (67) 
Al2o = 4L|2(9L3 - 9L2L1 + 2L?) (68) 
Aloi ~ 18Z/|Z/x — 18-ZjgZ/2 -|-27Z/2Z/2Z/^ — 13Z/2Z/^ — 4:SL^L'^Li — ^L^L/^L^ -\-\AL^L2L\ 

~2LsLl - 6Ll + lQ5LtLl/2 - 77LIl\ + SOL^L^ - QL2LI + L}72 (69) 
= -54L^ - 72LIL2L1 + A2LlL\ - ISL^Ll + UALaLjll - 7AL^L2L\ 

+22L3LV3 - QLlLi + IQLlLl - 20LlLl + IOL2LI - ALl/3 (70) 
Al2i = 81L2L2 - 27LlLl + IO8L3L2L1 - 9OL3L2L? + ISL3LI 

+271^/2 - 207 LlLj + ISGLIL^^ - 56L2L? + llLl/2 (71) 
Xlo2 = 36L|L2 + 36L^L? - 156L^L^Li + 52L^L2L? - 16L^Lf + 3L3L^ 

+93L3L^L? - 253L3L^Lt/3 + 35L3L2L? - 4L3L? - 41L^Li/2 

+135L^L?/2 - 62L^L^ + 76L^L73 - llL2L?/2 + Ll^ /2 (72) 
Aii2 = -144L^Li +48^21.2^2 +32L^L^ -30L3L^Li +126L3L^L^ -254L3L2Lf/3 

+10L3L^ +9L^/2 -43L^L2/2 +27L^L^ -25LIlI +67L2LI/6 -3Lf /2 (73) 
Ai22 = 216L^L2Li - 72LlLl - 27 L^Ll - A^L^LIlI - 9L3L2LI 

+9L3LI + 63L^Li - 22hLlLl + igiL^Lf - 59^2^^ + QL\ (74) 
Alo3 = 36L^ - 7SLIL2L1 + ^SLlL\ + 2ALIlI - 63LIlIlI - 20LIL2L\/2, - SL^L^ 

-47L3L^Li + 466L3L^L73 - 298L3L2Lf/3 + IOOL3L2L73 - llL3L?/3 

+Ll + 3L^L^/2 - ISL^Lf /2 - LIlI/3 + 7L^L73 - 7L2Lf/Q> + L^^/Q> (75) 
Ali3 = 18L^L2 - ISOL^L^ + 122LIL2L\ + 22L^L73 + 15L3L^ 

-86L3L^L^ + 92L3L^L^ - I88L3L2LV3 + 25L3L73 

-3L^Li + lAL^Ll - SLlLl - 2LIl[ + IIL2LV3 - 2L}V3 (76) 

Ai23 = -SIL^L^ + 270L^L2L^ - SlL^Lf + 108L3L^Li 

-282L3LILI + IOAL3L2LI - 22L3LI/3 - 911/2 + 33L^L2/2 

-53LlLl + GALjll - A9L2LI/2 + 17^^76 (77) 

Alo4 = 48L^Li +24L^L2 -174L^L2L^ +166L^L73 -91L2l3Li/2 +1141L2l2l76 
-189L^L2Lf/2 + 17L^Lj/2 + 5L3L74 + 27L^L\lI/2 - 3'37L^LIl{/Q 
+mL^LlLl/3 - I3L3L2LI/I2 - LsLiy2 + L^Li/8 - 17L^Lf/4 
+101L^L78 - 28L^Lj/3 + 61L^L724 - 5L2LiVl2 + Ll^/2A (78) 
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Ali4 = 6OLIL2L1 - 84LfL? + 33LlLl/2 - 24QLIlIlI/2 + TOSL^LaL^/e 
-l'dLlLl/2 - YJUL\Lxl2 + \mUL\L\l'i - ^2L^LIl\/?, 
-3IL3L2LI/3 + 5L3LI/2 - Ll/8 + 9LlLl/A - 67L^L^/8 
+14L^LV3 - LlLl/24 + 5^2^1712 - L^VS (79) 

Al24 = -hALlL2 + ISL^L^ + 54L^L^Li + 78L^L2L^ - 32L^L^ 

-27L3L^/4 + 27LsLlLl - 293L^LIL\/2 + I99L3L2L73 - 27L3L74 
+3L^Li/4 - 5L^Lf/2 + 33L^L72 - 2LIl[ - 3IL2LI/I2 + L]^ /2 (80) 

Alo5 = 6L^L2 + l2L\Ll - 19L^L^Li/2 - U^LIL2LI/3 + 29L^L72 + L^L^/2 
-bLlLlL\/A + 2467L2l^LV36 - 409L^L2L712 + 131L^L?/36 
-LgL^Li/e + 4L3L^L5 - 223L3L^L76 + 211L3L2l[/9 - I3L3L2L73 
+2L3L}V9 + i^2^?/8 - 11^2^f/6 + 69L^LV8 

-79L3lV12 + 137L2l}772 - L2LIV4 + L|772 (81) 
Ali5 = L^{-3QLlLi - 27LlLl/2 + 57LlL2Lj - IILIL\/2 + QL^LILi 

-I5L3LILI - 8L3L2L73 + L3LI - LlLl/2 + 3L^l72 - L2L73) (82) 
Al25 = -18L^ + 36LIL2L1 - ALlL\ - 3LlLl/A - UlLjLlLl/A 
+111L^L2L74 - 239L^L736 - 31^1^1^ + I93L3LILI/6 
-87L3LILI/2 + 33L3L2LI/2 - 3IL3L7I8 - L78 + 2LlLl 
-81LiLi/8 + 163L^Lf /12 - 45L^L78 + 3L2Li74 - L}772 (83) 

where the invariant traces Li, L2 and L3 are defined as in the main text Eq. 15. Some 
useful supplementary polynomials (used here and in the main text Eqs. 47-52) are: 

Ls := (Li - L3)/3 = (me + m^)(m^ + m^)(m^ + me) (84) 
Lp2 := (3L2 — Li)/2 = mg + m^ + m^ — mgm^ — m^m^- — m^-mg (85) 
Lp3 (3L3-L2Li)/2 = m3 + m3 +m3 

— {mlm^ + m^me)/2 — {m^^rrir + m^m^)/2 — {m^me + m1mr)/2 (86) 
Lp6 Lp^ — L^/A — (mg + + m^ — m^rrin — rrii^mr — m^nie) 

X (mg + + - mgm^ - m^m^ - rriT-ml) (87) 

where the charged-lepton discriminant L\ is defined in the main text, Eq. 17. 
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